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ABSTRACT 

Observations of the CO III] 88.35um line and the [01] 
63.2um line have been made with a new far- infrared spectrom- 
eter. The sources M17, NGC 7538 and W51 have been mapped in 
the [0 III] line with 1 arcminute resolution and the emission 
is found to be quite widespread. In all cases we find that 
the peaK of the emission coincides with the maximum radio 
continuum. Where possible we have simultaneously mapped the 
far- infrared continuum and in M17, NGC 7538 and W51 the 
continuum peak is found to be distinct from the center of 
ionization. The [0 III] line was also detected in W3 , W49 
and in a number of positions in the Orion nebula. Upper 
limits were obtained on NGC 7027, NGC 6572, DR21, G29. 9-0.0 
and M82. The 63.2pm [01] line has been detected in N17, 

M42, and marginally in DR21. A partial map of M42 in this 
line shows that most of the emission observed arises from 
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the Trapezium and from the bright optical bar to the 
southeast. The transition of HD at 112ym was also 

searched for in the direction of the Kleinmann-Low nebula 
but without success . 


I. INTRODUCTION 

Far-infrfikred fine structure lines provide a peftieuldt^ly 
useful means of studying ionic species in H II regions* plane- 
tary nebulae and galaxies. Development of this field has only 
become possible in recent years with the provision of high 
altitude telescopes such as the Kuiper Airborne Observatory. 

The first detection of [0 III] X = 88.35um was by Ward et al. 
(1975) in M17. Since then, the same line has been detected 
in Orion (Baluteau et al. 1976 and Moorwood et al. 1978) and 
in several other H II regions (Dain et al. 1978). So far the 
sources have either been observed with large beams, or the 
signal-to-noise ratio has been too low to produce detailed maps. 
We present in th'.s paper the first far-infrared fine structure 
line maps, which were taken at a resolution of 1 arc minute. 

In general, the distribution of [0 III] emission can be 
expected to resemble that of the radio continuum (since both 
arise from electron collisions) but with some important dif- 
ferences. First, the radio concinuum (free-free emission) is 
proportional to the emission measure or electron density 
squared, while the fine structure line is, except at very low 
densities, proportional to ionic column density and hence in 
most cases varies linearly with the electron density. Sec- 
ondly, the abundance of O'*"*’ relative to hydrogen is not nec- 
essarily constant from cloud to cloud, or indeed within a 
cloud. Apart from variations in the local oxygen abundance, 
molecules and grain surfaces can lower the effective atomic 
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oxygen content of the ^as._ Furthermore, the intensity 
emd effective temperature of the ionizing radiation will 
affect the relative abundance of 0^^. Because 3S.11 eV is 
required to ionize oxygen to 0^^, and SS.89 eV to ionize it 
further, 0^^ is an abundant species in many H II regions 
but may sometimes appear depleted when the exciting star is 
cool. Comparison of radio and fine structure maps, partic- 
ularly where differences exist, thus provides information on 
the spectral type of the stars underlying the H II region, 
and may possibly shed some light on their formation. Com- 
parison with fine structure lines of other ions would also 
yield important information- An interesting example is 
[S IV] at 10. 5y. Although this ion has a very similar ioniza 
tion requirement (34.83 eV) to [0 III], much higher electron 
densities are needed to reach thermal equilibrium. 

The [0 I] line at 63.2ym is of particular interest 
because it can be excited not only by the ionized material in 
H II regions, but also by warm, neutral atoms in post-shock 
gas (Hill and Hollenbach, 1978). The CO I] line has recently 
been detected in M42 and M17 by Melnick et al. (lS79a) but 
the large beam-size used (4x6 arcminutes) leaves unanswered 
the question of what kind of conditions are producing the 
emission. We have confirmed the presence of [0 I] in M17, 
and in M42 examined each of the regions most likely to 
produce the emission. 

Detection of HD in molecular clouds would be a very 
direct way of determining the D/H ratio in these regions. 
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Several deuterated molecules have already been detected by 
radio techniques, but uncertainties about the exact nature 
of the chemical fractioxtation processes involved in the for- 
mation of these molecules make the results extremely diffi- 
cult to interpret. However, because molecular hydrogen is 
the principal species in these clouds by such an enormous 
margin, the HD/Hj ratio is an accurate representation of the 
cloud D/H ratio. In addition, rotational transitions of HD 
may be significant in the overall cooling of molecular clouds. 
A preliminary search has hence been undertaken for this 
molecule . 


II. INSTRUMENTATION AND OBSERVATIONS 
The spectrometer uses a fixed wavelength Fabry-Perot 
interferometer cooled to 4.2 K in conjunction with a high 
order scanning Fabry-Perot at room temperature. The scanning 
Fabry-Perot, although warm, is placed immediately outside the 
cold aperture stop in such a way that at wavelengths outside 
its narrow passbands it reflects only liquid helium ■ emperatures 
into the detector. High resolution is thus obtained while 
minimizing the thermal background radiation on the detector. 
Because only a narrow range of wavelengths is incident on the 
detector at any one time, "guiding" noise and fluctuations in 
the atmospheric water vapor content along the line of sight 
do not contribute significantly to the system noise. The 
Fabry-Perot mirrors consist of gold-plated nickel mesh 
stretched tightly over invar rings; each filter has a finesse 
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of mbout 3 S. Short -wavelength radiation is blocked by a 
2 lam thick BaF 2 filter (AR coated) and a quartz/dieu&ond 
8 catter/CaF 2 filter, both at 4.2 K. The detector is a 4.2 K 
gallium-doped germanium photoconductor mounted in an integrat- 
ing cavity, both kindly lent to us by K. Shivanandan. 

The spectrometer resolution was determined in the 
laboratory by scanning across the output of a far-infrared 
methanol laser operating at 96. Sum and at 118. Sum. The 
response profile is a lorentzian (as expected for a Fabry-Perot 
system) with a full width at half maximum of 0 . 12 um when the 
system is used in 20 th order. 

The observations were made with the 91.4 cm telescope of 
the Kuiper Airborne Observatory in four flights during 
September, 1978, at an altitude of 41,000 feet (12.5 km). 

An oscillating secondary mirror chopped the beam in azimuth 
by between two and six arcminutes (depending on the source) 
at a frequency of 35 Hz. A typical observation consists 
of a 40 second scan with the source in one beam followed by 
an equal time with the source in the other beam. The 
detector output is amplified by a room temperature JFET 
amplifier, then synchronously detected at 35 Hz. After 
passing through a six-pole low pass Bessel filter the out- 
put is sampled 1024 times each scan and the result stored 
on magnetic tape for later analysis. 

Wavelength calibration was performed in flight by 

18 

observing a gas cell containing HjO, H2 0 or HDO, chopped 
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ag ains t a. IQQ^C hlaekbody — bl«c)cbody was monitored 
frequently during observations to check for drifts in sensi- 
tivity. Absolute intensities vrare calibrated by observation 
of Jupiter and by the expected continuum of the nebulae them- 
selves. The absolute calibration uncertainty is of the order 
of 30%, although when mapping a given source the relative 
uncertainty across the map is much less. The measured beam- 
width was 1 arcminute (full width at half maximum) but obser- 
vation of extended sources shows that the effective area’ of 
the be£un was equivalent to that of a 75 arcsecond disc. 
Pointing was determined by offset from guide stars and is 
accurate to within 15 arcseconds. 

At 88um the instrumental resolving power (X/6X) was 
700, with a system NEP in flight of 1.4 x 10 ”^^W//Hz (includ- 
ing all telescope and spectrometer losses). At 63 and 112ym 
the performance was not as good, due partly to mistuning of 
the 63um filter, and partly because the scanning Fabry-Perot 
was operated in lower orders, limiting the resolution. 

Ill . THEO.RY 
a) [0 III] X = 88.35um 

Russell-Saunders coupling is a good approximation for 
the ground term of 0 , which is " ^210* Levels 

of different J in this term are split by the spin-orbit inter- 
action into a normal triplet (see Figure 1). Since all of 
these levels have the same parity, there can be no electric 

dipole transitions between them. The lowest-order transitions 

3 3 

are the magnetic dipole (Ml) transitions ?2 and 
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3 3 

Pq, and the electric quadrupole (E2) transitions 
3 3 3 3 

?2 Pq end P 2 Pj^. The Ml transitions, having Einstein 

-5 -1 

A-coefficients of “vlO sec , are the brighter by factors of 

about 10^. These transitions are easily excited and de- 

excited by the magnetic fields of passing electrons. 

Typical electron densities of H II regions are in the 
-3 4 -3 

range 400 cm 5 S 10 cm , with temperatures in the range 
3 4 

5 X 10 K s T < 10 K. Under such conditions, the populations 
of the fine structure levels of the ground term of 0 ** are 
close to the thermal equilibrium values, which in turn are 

close to the statistical weights for the states because 

kT >> E. The radiation from the ionized clouds in the 

3 3 3 3 

?2 - 51.8vm) and P^ ^q * 88.35um) transitions 

is optically thin. These facts imply a very simple expres- 
sion for the power per unit area and solid angle in these 
lines: 


^(k>j) ^ 471 ^ 

where is the A-coefficient for the transition, n^^ 
the density of o"*"*’ ions in the upper state k, and I is the 
length of the nebula along the line of sight, n^^ can be 
related to the total density n of O'*"*’ ions by 


where k = 1,2, and f^ 
3 

ratios of the and 


n 


k 




. ^1 ^ - '^2 
- — and fr, = — 

n n ‘ n n 


are the population 


?2 levels respectively, relative to 
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the population of Pq. We are ignoiing the populations of the 
more highly excited electronic states, which are negligible 
at these temperatures and densities. By using the principle 
of detailed balance and knowing the excitation cross-sections 
for e" - 0^^ collisions (Spitzer, 1978), the quantities f^^ 
and f 2 can be calculated as functions of electron density 
and temperature, as shown in Figure 1. Thus, Equation (1) 
has two uses; (a) with an assumed 0^"*^ abundance and a value 
for n^, one can predic.t the strengths of the 51.8ym and 
88.35ym lines; or (b) with a measurement of the strength of 
either line one obtains directly a value for the O"*”'’ column 
density. Inspection of Figure 1 shows that the population of 
the level as a fraction- of the total population is relatively 
constant over a large range of electron densities, and hence 
observation of the 88.35ura transition alone gives an accurate 
measure of the O’*"*’ column density, even when large uncertain- 
ties exist in the assumed electron density. Relative abun- 
dances of two ionic species which both have far-infrared fine 
structure lines can hence also be well determined by measure- 
ments of this type. 

To derive a fractional abundance of 0^ ' relative to 

hydrogen, an accurate value for the electron column density 

is required. The derivation of electron column densities 

from radio continuum data is difficult, however, because it 

is necessary to model the source with a shape and structure 

in order to determine the relation between electron column 

2 

density (n^l) and emission measure (n^ 1). Only in very 
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simple cases can the electron density be determined 
accurately from the emission measure, and this represents a 
major source of uncertainty when fractional abundances are 
calculated on this basis. 

An alternative approach is to measure the intensities 
of both the CO III] lines, allowing the electron density to 
be determined in a more direct manner. From detailed balance 

3 

in the P term, the following relation can be easily derived; 
n = ^01(^21 ^20) * ^02(^21 ~ a ^10} (2) 

^0l(i ^12 " ^21 " ^20) * ^02(1(^12 * ^10) " ^21) 

Inhere y^j is the collisional rate coefficient for the transi- 
tion i-^j (see Spitzer 1978), A^j is the A-coef ficient for 
the transition i-»-j , a = ( 51. 8/88 . 35 ) (Aj^q/A 2 j_) = 0.158, and 
F = I(88.35ym)/I(51.8ym) is the ratio of the observed line 
intensities. A procedure similar to this for determining 
has already been employed by Melnick, Gull, and Harwit (1979b) 
and produces values of n^ significantly different from the 
radio-determined values. While this method is more direct, 
it suffers from two difficulties. First, the ratio of intensi 
ties, R, does not vary linearly with n^ but for a given frac- 
tional change in density changes very slowly at both the high 
and low ends of the range of typical densities. Thus the 
derived value for n^, which represents an average of the 
dependence upon n^ .'•f R over the line of sight, should be a 
lower limit to the average value of n^ and may be considerably 
lower than the true value. It is thus anomalous that some 
of the values for n^ derived via this method by Melnick et al. 
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are higher than the rms values of obtained from emission 
measures, which should be upper limits to the average value 
of n . Secondly, conversion of n. to a column density requires 
knowledge of the line of sight length of the emitting region, 
and hence the size of the object must be well known. 

b) [0 I] X = 63.2ym 

2S+1 3 

The ground term of neutral oxygen has " ^2 1 0 ’ 

just like 0^^, but the fine structure triplet is inverted as 
shown in Figure 2. As before, the most important transi- 

0 0 it 

tions are the Ml transitions, Pq -► (145. 4um ) and 

3 3 

Pi -► p 2 (63.20vm). The collisional excitation processes 
for 0 I are very different from 0 III, though, for t;.o impor- 
tant reasons. First, 0 I is neutral, so the cross-sections 
for electron collisions will not be as large and electron 
collisions will not be as predominant an excitation mechanism. 
Secondly, 0 I has an ionization potential of 13.618 eV compared 
to 13.598 for hydrogen, so much of the neutral oxygen will be 
outside of the H II region and depend upon collisions with 
hydrogen atoms for excitation. The fact that 0 I is a con- 
stituent of the H I - H II transition region makes these fine 
structure lines very important for studying the structure of 
ionization fronts. It also makes prediction of the strengths 
of these lines very difficult. It is beyond the scope of 
this paper to study this problem in detail*, we can, however, 
make a few reasonable estimates. 

*taken from the data presented in Bashkin and Stoner (1976), 
which is more accurate than the compilation of Moore (1949) 
that led to the previously-accepted value of 147um. 
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Ex«alnatton of the colllsienal rate ooefficients for 
oxygen oollieione with protone (Batwaall and Wolf, 1936), 
electrone (Saraph, 1973) and hydrogen atone (Launay and 
Roueff, 1977) finde the electron ratee to be enaller than, 
but competing with, the proton ratee. The hydrogen rates 
are down a factor of about 10 from the other two for the 
same temperature and density. Since the hydrogen collisions 
will contribute only in the lower-temperature non-ionized 
region, it is possible that the emission from this source 
will be a factor of 'v«100 down from the emission from the 
ionized side. Taking electrons and protons to dominate the 
excitation, we can produce the population ratios in the P 
term shown in Figure 2. Prediction of the [01] intensities 
then follows directly from an assumed oxygen column density. 

Ill . RESULTS 
a) tO III] 

We were able to observe only the 88.35um line in this 
flight series. Figure 3 shews the result of a pair of obser- 
vations; each represents two 40-second scans of the indicated 
wavelength range, or a total of about 15 seconds of integra- 
tion on the line itself. The profiles were fitted by a 
least- squares procedure with a lorentzian line shape, which 
is a good approximation to our instrumental response. We 
therefore report peak intensity rather than integrated flux, 
as it is assumed that the line is substantially narrower than 
our resolution. At 38.35um our resolution was sufficient to 
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bring- the Co IIIl line well clear of all atmospherle' water 
vapor lines, and no correction needs to be made for absorp- 
tion. 

Because the fine structure line and the 8 8 urn continuum 
are recorded simultaneously by the spectrometer, our maps 
have essentially perfect registration between these two 
components. One expects the peak of the fine structure 
emission to be coincident with the maximum radio continuum 
emission, and this appears to be the case in all the sources 
we observed. We can thus determine the position of the far- 
infrared continuiom peak in these sources with some certainty, 
relating it to radio peak positions which are usually known 
to high accuracy. 

Table 1 is a summary of our intensities for the brightest 
peak in each source, together with predicted intensities and 
derived O'*"*’ column densities. In the latter two quantities 
we have used an electron temperature of 8000 K, with electron 
densities and sizes of H II regions from the 15 GHz, 2 arc- 
minute resolution work of Schreanl and Mezger (1969) — except 
for G2 9. 9-0.0, for which we referred to Zeilik (1977). The 
predicted intensities for che H II regions assume an O'*’"*’ 

_4 

abundance of 1.5 x 10 relative to hydrogen. This is 

. . -4 

derived from a cosmic oxygen abundance of 6.5 x 10 , of 

which 2 3% is assumed to be in the form of The predic- 

tions for the observed planetary nebulae NGC 7027 and 'IGC 6572 
are due to Simpson (1975). 
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In principle it is possible to use our 0^^ column 
•densities to derive values for the 0** fractional abundance 
by comparing them to the electron column densities. However, 
electron densities are subject to a number of uncertainties, 
as discussed in Section XIa, so that reliable values of the 
0 fractional abundances cannot yet be derived. In Table 2 
we show first the electron densities and linear sizes of the 
nebulae derived from radio data, and the 0^^ fractional 
abundances derived from these estimates. In the next three 
columns are the electron densities from Melnick et al. (1979b), 
the nebular sizes from this present work, and the 0 frac- 
tional abundances derived on this basis. In addition, the 
fractional abundances obtained by Melnick et al. are listed. 

It is clear that substantial differences result from the use 
of different techniques for estimating electron column densi- 
ties. This problem is likely to remain until more detailed 
models of the nebulae can be produced. It is also apparent 
that our values for the 0^'*’ abundance are considerably lower 
than those of Melnick et al., who used source dimensions based 
on radio data and hence underestimated the size of the emit- 
ting region. We thus do not find it necessary to argue for 
an increase in the commonly accepted, oxygen abundance ratio. 

Figures ‘•^-6 are the contour maps we have obtained. In 
each map, the positions where data was taken are marked by 
crosses; the data at each point was usually in the form of a 
pair of 40-second scans. Solid contours show the [0 III] 

\ = 38,35um line intensity in units of 10” cm”'’, and the 

dotted contours show the continuum intensity in our 0.19um 
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effective bandwidth (ir/2 times the resolution), in the same 
units . 

In each case studied we find that the far- infrared 
continuum is not centered on the peak of the ionization. 
Although this is not an unexpected result, it suggests that 
simple spherical models of H xl regions which have the dust 
uniformly mixed with the ionized material will not give 
realistic predictions of the detailed properties of these 
sources. 

In most sources examined the fine structure line intensity 
is distributed rather widely and smoothly. This is in part 
because our beam width of 1 arcminute is larger than some of 
the fine spatial structure, but in part must represent the 
actual widespread and relatively smooth distribution of 
[0 III ’ in regions where the density is rather low and the 
Stromgjen spheres of the ionizing stars are very large and 
tend to overlap. 

•v'e report below the details of the results from individual 
objec js. 


i) M17 

In This source we mapped the far-infrared continuum 
simultaneously with the 88.35ym line, as shown in Figure 4. 

The peak [0 III] emission is centered at the peak of the radio 
continuum (G15.0-0.7) and closely follows the radio contours, 
exter.alng for several arcminutes along a northwest- southeast 
line." le also looked for and detected CO III] at one Dosition 
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in the bright optical region to the northeast, at an intensity 
-16 —2 

of 1.9 X 10 W cm . This position is close to tha peak B 

observed by Montgomery et al. (1971) at 3.4 tm. No emission 
-17 -2 

(<2x10 Wcm ,3o) was seen from the 90 pm northern 
source of Lada (1976). The peak [0 III] intensity is in 
reasonable agreement with the value predicted from the radio 
data. The total [0 III] intensity is also consistent with 
that measured by Dain et al. (1978) in a 4* x 4.4* beam. 

It is clear that the CO III] emission is very widespread, 
and probably extends well beyond the limits of our map. This 
supports the contention of Lada (1976) that a large group of 
stars have formed along the edge of the molecular cloud and 
are supplying an enormous flux cf ionizing photons. 

The 88um continuum radiation forms an extended ridge 
along the molecular cloud to the southwest of the [0 III] 
peak. This is similar to the recent findings of Wilson et al. 
(1979) who suggest that thermal emission from dust at the 
interface of the ionized region and the molecular cloud is 
responsible. We find that the continuum emission peaks at 
the northern edge of the ridge. Wilson et al. do not see 
such a peak at S9um, and it is possible, although unlikely, 
that the discrepancy is due to the different wavelength used. 
An earlier map by Harper et al. (1976) has been made at lOOum, 
but at their resolution of 2.3 arcminutes the continuum peak 
appeared coincident with the radio source GlS.0-0.7. 
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- - - - - - ii) N6C 7538 

This is the first detection of the [0 III] line in this 

source. The observed peak intensity is in good agreement with 

the value predicted from radio studies, assuming an 0^* 

•4 

fractional abundance of 1.5 x 10 . 21 cm interferometry by 

Israel et al. (1973) shows that the ionization resembles a 
section of a spherical shell centered on a point to the 
northeast of the peak emission. The CO III] peaks in the 
same place as the radio continuiim and shows a similar struc- 
ture (see Figure 5) though we do not resolve the bright con- 
densations. We observe maximum far-infrared continuum at 
the position of IRS2 (Wynn-Williams et al. 1974) , in agree- 
ment with recent work by Werner et al. (1979). At this 

—17 -2 

position no CO III] is observed (<2 x 10 W cm , 3o) , 

-17 -2 

consistent with the upper limit of 4 x 10 W cm of 
Moorwood et al. (1978). 

iii) W51 

The results of observations in this source are shown in 
Figure 6. As expected, the maximum CO III] emission comes 
from the region of most intense 5 GHz radio continuum, desig- 
nated G49.5e by Martin (1972). This position is coincident 
with the infrared source IRSl (Wynn-Williams et al. 1974). 

The far-infrared continuum peaks somewhat to the north and 
west of this source, close to the maximum 55ym continuum 
observed by Harvey et al. (1975). In contrast with the SSym 
result, however, we find the far-infrared continuum to be 


quite extended, implying a widespread distribution of cooler 
dust throughout the ionized material. Although we cannot 
rule out the possibility that the far-infrared peak is centered 
on IRS2, it does appear more likely that it lies between IRSl 
and IRS 2. 

Using the electron density and linear size that Martin 

—16 —2 

derives for G49.5e, ve can predict a value of 0.9 x 10 W cm 
in the CO III] line for this peak, in better agreement with 
the observations than the prediction in Table 1 , based on the 
lower resolution 15 GHz data. 

iv) W49A 

Although we were able to detect [0 III] in this source, 
no attempt was made to map the spatial distribution. Refer- 
ing to Table 1, it is clear that cur observed value for the 
peak flux is considerably lower than that predicted from the 
radio work of Schraml and Mezger. This may be because much 
of the flux seen by Schraml and Mezger is due to bright con- 
densations (Wynn-Williams, 1971 ) too small for us to resolve. 

2 2 

This would imply that <n^ > > 60(<ng>) , i.e. rather extreme 
conditions with ionization condensations filling perhaps 1Q% 
of the beam. Another possibility is that the underlying 
stars are too cool to ionize much of the oxygen to and 

hence that the fractional abundance of this ion is low. It 
would be interesting to search for the 10.52um line of CS IV] 
(required ionizing energy 34.8 eV; essentially the same as 
for 0 III) to test this possibility. 


W3A is unresolved by our beam; the peak CO III] emission 
seems to be in the direction of W3.-IRS1 (Wynn-Williams et al. 
1972). Negligible [0 III] was observed in positions 1* north, 
1’ south, 1' east and 1* west of the peak, or from IRS4. 

Our observations of IRS3 show no emission, but the upper 
limit is only about half that of the detection in IRSl. 
Unfortunately IRS3 was only observed at the beginning of the 
flight series, when the spectrometer was adjusted for low 
resolution and hence unable to discriminate well against the 

continuum of the source. The total intensity we see 

—16 —2 

(1.0 X 10 W cm ) is considerably lower than that seen by 
Dain et al. (1978) (3.1 x 10"^®W cm"^), but higher than the 
upper limit of Moorwood et al. (1978) (<7 x lO'^^W cm ^). 
Although W3A appears to be a point source, it is possible 
that the large beam used by Dain et al. included emission 
from regions we did not examine, "his seems unlikely as all 
such regions have lower radio continuum flux than IRSl. 

vi) M42 (Orion) 

A number of positions in this source were observed, as 
shown in Figure 7. The emission peaks in the vicinity of 
the Trapezium, with an intensity of 1.2 x 10 cm ^sr 

in good agreement with the value of 1.1 x 10”^'^W cm”^sr ^ 
observed by Moorwood et al. (1978). These authors found 
that the decrease in emission intensity with increasing dis- 
tance from the Trapezium was very gradual. The present data, 
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at & ipdsittSHS preseni: further 

" ■ for distribution of tO III] emission, 

see SSI' intensity of 70% of the peeik value at eadh of the 
four poiistir 2 arcidhutes away from the Trapesium, and a con- 
tinued ^adual fall-off at tn4 pos it ions to the northeast 
and sWth^st ;''- ' 

" Altogethef' tlie distributi^ eiaission is considerably 
more xuilfdihn "fbAn that pt^dicted by t^ 

C1973). T^ire it also a significant difference between our 
data fpoin^ 2 arc minutes north and east and 2 arc minutes 
south 'and west, showing that even at 1 arcminute resolution 
the assumption of spherical symmetry for the source is not 
very realistic. 

Under the assumption of a homogeneous cloud with uniform 
ionization fraction throughout, the CO III] line intensity 
should be approximately proportional to the square root of 
the radio continuum flux. Comparison with the 15 GHz data 
of Schraml and Mezger (1969) does show some such correlation, 
but also significant differences. Variations in the density 
of the ionized gas are most likely responsible for these 
differences. Qualitatively they correspond to the presence 

_3 

of low density (n^ > 100 cm ) material along the line of 
sight, presumably in front of the denser ionized region which 
produces most of the radio continuum. Since M42 is so large 
and time on the source was limited, a complete map on the 
scale of our beam size was not attempted. It is clear that 



21 


a ditilleid the Orion nebula in the CO IIIl line 

would provide a valuable probe of low density iohi^ed i^terial 
and of ionization conditions, particularly in the obscured 
regions such as the dark bay to the east. 

b) CO 13 

Observations of the [01] line were hampered by a Lower 
system performance at this wavelength. The flight schedule did 
not allow observation of Jupiter at 6 3 urn, and hence for cali- 
bration we must rely on observations of the moon and of the far- 
infrared continuum of the sources themselves. In addition, there 
is a strong atmospheric water vapor line at 63. 3 urn. Although 
this line is quite narrow and does not attenuate the [01] 
emission, it absorbs an unknown amount of the continuiun radi- 
ation adjacent to the [01] line and hence produces an appar- 
ent decrease in the line intensity. This effect is most 
important for the stronger continuum sources . The atmospheric 
water vapor content can change by a factor of 3 or more during 
a flight, and is a strong function of altitude. Because 
of these uncertainties our estimates of the absolute [0 I] 
intensities are limited to an accuracy of about 50%. 

Ionospheric [01] emission at 63.2ym has recently been 
detected by Grossman and Offermann (1978). Although the 
intensities they see are comparable to our observed values , 
our use of beam switching throughout the flights eliminates 
any possibility of contamination from this source. In 
principle an optically deep ionospheric line could severely 
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_ hpMeveag ,o alc t il* t iO B r- •hows - 

-;,^-^ii^r^.^to -b«. pr^blaa. 

; ^ Ill Table i the observed CO 1 3 intensities t 

together with predicted intensities based on the radio data 
of Schraml and Mesger (1969) and an assuned fractional 
oxygen abundance of 1.5 x 10 . The predicted CO I3 intensi- 

ties are in qualitative agreement with our results, although 
too much reliance should not be placed on the predictions. 

The excitation of neutral oxygen in H I - H II transition 
regions is a complex problem that deserves a great deal more 
theoretical attention, as well as additional observations. 

A detailed discussion of CO I] emission from weak D-type 
ionization fronts and shocked molecular clouds is available 
in a recent paper by Hill and Hollenbach (1978). 

i) M17 

We have detected the line in the direction of the peak 

-16 -2 

radio emission in M17, with an intensity of 1.2 x 10” W cm 
Although no attempt was made to map this source, comparison 
with the result of Melnick et al. (1979a) of 24 x 10“^®W cm“^ 
in a 4 X 6 arcminute beam (an area 15 times larger than ours) 
suggests that the emission is widely spread throughout the 
nebula. 

ii) DR21 

This is the first detecxion of a far-infrared fine 

structure line in this source. The observed flux is 
-17 -2 

8.1 X 10 W cm , and although random noise is low, the 
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difficulties of correction properly for the atmospheric vapor 
prevent us from ascribing more than a 3c certainty to this 
detection. 

iii) M42 (Orion) 

^ throughout the nebula were observed, as 

shown in Figure 8. It is apparent that the majority of the 

[01] emission is from the ionized regions, or more probably 

from the ionization fronts associated with the H II regions. 

This is shown in the strong emission coming from the bright 

optical bar to the southeast, which also shows up strongly in 

the 6300 A line of [0 I] (T. R. Gull, 1978). The strongest 

“16 “2 

signal we see is at the Trapezium, 1.3 x 10 W cm 
Assuming that the emitting region uniformly fills the effec- 
tive 75 arcsecond disc of our beam, this corresponds to a 

-1 -2 

flux at the surface of the nebula of 0.039 erg sec cm 
This is comparable to the values predicted by Hill and 
Hollenbach (1978) for the total ionization front plus shock 
front emission following a 5-10 km s ^ shock. However, since 
no molecular hydrogen radiation is seen from regions strong 
in Co I] emission nor is much CO I] emission seen from around 
KL where the vibrationally excited H2 emission is most intense, 
it seems clear that shocks are not the main mechanism for 
exciting CO I]. The expected alternative location for form- 
ing the Co I] line is in the transition zone between H II 
regions and the neutral material; these are in fact the regions 
we observe to be most intense. 
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If the Trapezium is indeed the brightest source of 
Co I] emission, then our observed surface brightness is in 
conflict with the value of 0.13 erg sec** cm’ , averaged over 
the whole nebula, determined by Melnick et al. (1979a). Two 
possible explanations are (a) the majority of the emission 
comes from regions we did not observe, or Cb) our chopper 
throw of 5 arcminutes was insufficient to bring our reference 
beam clear of the emission. Neither explanation appears 
attractive, however, and it is more likely t^at the difference 
is due to the great difficulty of correcting properly for 
atmospheric absorption so close to a water vapor line. 
Observations at higher spectral resolution than has so far 
been used would resolve this difficulty in a straightforward 
manner. Even allowing a factor of two error in each of the 
two measurements, a very widespread distribution of [0 I] is 
implied by the large amount of flux in the 6x4 arcminute 
beam relative to the 75 arcsecond beam. This is further 
evidence against a shock mechanism being the primary source 
of excitation of the line. 


c) HD 

_4 

The HD molecule has a dipole moment of 5.9 x 10 D 
(Trefler and Gush, 1968), resulting in a pure rotational 
spectrum, the J=l-^0 transition of which falls at 112.08um. 
The Einstein A-coeff icient of this transition is 2.54 x 10’ 
sec’^, and the transition is thus in thermal equilibrium at 
quite low densities. Calculation of expected line intensi- 
ties is similar to the case of atomic fine structure lines, 
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principal collision 

partner. For exanplet Bussoletti and Stasinska (1975) have 
predicted the flux from the Orion molecular cloud as 

I. 1 X 1(J“^®W cm“* for,the 112pm l^iie. 

We searcifed for this transition in the direction of the 

KleiiraMinn-LoW nebula , but after five minutes of integration 

the baseline stability became the limiting factor, due to a 

combination of the lower resolution used at this wavelength 

CX/dX * 180) and the high continuum level of the source. The 

—17 -2 

upper limit we obtain of x 10 W cm in our beam is not 
sufficiently low to rule out any but very optimistic models 
of the Orion molecular cloud. On the other hand, increased 
spectral resolution and other feasible improvements to the 
spectrometer should allow detection of this and other far- 
infrared molecular lines in warm clouds. 
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Table 3 

Suonary of [O l] Resxilts 


Source 

HI 

£(1950) 

Observed Co ll 
Xatenslty 

c*”^) 

Predicted CO I] 
Intensity 

(10“^Sl c*“^) 

M17 

18**17^f5 

-16®13’24" 

1.2 

2.0 

lf42 

5**32*49f6 

-5®25*16" 

1.3 

1.9 

DR21 

20**37*13?5 

+42®C3*5l'* 

0.81 

J 

0.60 
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FIGURE CAPTIONS 

Fig. 1. Relative populations of the levels of the ground 

term of 0^* plotted as a function of electron density 
for four different electron temperatures. The 
inset shows the energy level diagram for the ground 
term of this ion. Comparison of this result with 
that of Smith (1969) shows that the effect of the 
higher-energy electronic states is negligible in 
this density range. 

Fig. 2. Relative populations of the levels of the ground 
term of atomic oxygen plotted as a function of 
electron de'ssity for three different electron tem- 
peratures. The inset shows the energy level diagram. 

Fig. 3.. [0 III] emission in M17 and NGC 7538. Each profile 

represents one pair of 40 second scans. A lorentzian 
profile has been fitted to the data via a least- 
squares procedure, yielding values for both the peak 
line intensity and the continuum level. 

Fig. 4. Map of [0 III] emission in M17 (solid contours) and 

88um continuum (dotted contours). The h marks a 

position on the bright optical bar where [0 III] was 

— X6 *2 

detected at an intensity of 1.9 x 10 W cm 

Fig. 5. Map of [0 III] emission in NGC 7533. The A marks 
the position of IRS2 (Wynn -Williams et al. 1974 ). 
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Fig. 6. Map of Co III] emission in W51 (solid contours) 

and 88um continuum (dotted contours). The positions 
of IRSl and IRS 2 are from Wynn-Willicuns et al. (1974) 
Fig. 7. The Orion nebula, showing positions at which CO III] 
data was taken with measured intensities in units of 
lO’^^W cm“^. 

The Orion nebula, showing CO I] data with intensities 
in units of 10”^^W cm”^. 


Fig. 8. 
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